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We study noise and transport in multiterminal diffusive conductors. Using a Boltzmann-Langevin 
equation approach we reduce the calculation of shot-noise correlators to the solution of diffusion 
equations. Within this approach we prove the universality of shot noise in multiterminal diffusive 
conductors of arbitrary shape and dimension for purely elastic scattering as well as for hot electrons. 
We show that shot noise in multiterminal conductors is a non-local quantity and that exchange effects 
can occur in the absence of quantum phase coherence even at zero electron temperature. It is also 
shown that the exchange effect measured in one contact is always negative - in agreement with 
the Pauli principle. We discuss a new phenomenon in which current noise is induced by thermal 
transport. We propose a possible experiment to measure locally the effective noise temperature. 
Concrete numbers for shot noise are given that can be tested experimentally. 
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I. INTRODUCTION 



Shot noise is a noncquilibrium fluctuation of the cur- 
rent in mesoscopic conductors caused by random flow 
of the charge. It can be thought of as an uncorrelated 
Poisson procesaj giving rise to a simple formula for the 
spectral density of the shot noise, S"^ — el, where / 
is the current through the conductor and e is the elec- 
tron charge. Being the result of charge quantization, the 
shot noise is an interesting and highly nontrivial physical 
phenomenon. □ In contrast to the thermal fluctuations of 
the current, the shot noise provides important informa- 
tion about microscopic transport properties of the con- 
ductors beyond the linear response coefficients such as 
the conductance. For instance, the shot noise serves as a 
sensitive tool to study correlations in conductors: while 
shot noise assumes the Poissonian value in the absence of 
correlations, it becomes suppressed when corifelations set 
in as e.g. imposed by the Pauli principle.0Ll In partic- 
ular, the shot noise is completely suppressed in ballistic 
conductors ,□ and it appears thus only in the presence of 
a disorder. 

In diffusive mesoscopic two-terminal conductors where 
the inelastic scattering lengths exceed the system size the 
shot-noise suppression factor for "cold" electronsii.e. for 
vanishing electron temperature) was predictedHliJ to be 
1/3. The suppression of shot noise in_diffiusive conduc- 
tors is now experimentally confirmed. E3~t£l While soma 
derivations are based on a scattering matrix.ai>i)roachB'tll 
or conventional Green's function techniquetda and thus 
a priori include quantum phase coherence, no such effects 
are included in the semiclassical Boltzmann-Langevin 
equatipaj-approach, which nevertheless leads to the same 
result £3'Ej However, while in the quantum approach for 
a two-terminal conductor the factor 1 /3 was even shown 
to be universal,li3 the semiclassical derivations given so 
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FIG. 1. a) Multiterminal diffusive conductor of arbitrary 
2D or 3D shape and with arbitrary impurity distribution. 
There are A*' leads with metallic contacts of area L„, and 
Im is the m-th outgoing current. S denotes the remaining 
surface of the conductor where no current can pass through, 
b) Conductor of a star geometry with N long leads which 
join each other at a small crossing region. The resistance of 
this region is assumed to be much smaller than the resistance 
of the leads, c) Wide conductor: the contacts are connected 
through a wide region, so that the resistance of the conductor 
comes mainly from the regions near the contacts, while the 
resistance of the wide region is negligible, d) H-shaped con- 
ductor with four leads of equal conductances, G/4, connected 
by a wire in the middle of conductance Go. 
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are restricted to quasi-onedimensional conduc- 
tors. Thus, although phase coherence is believed, not 
to be essential for the suppression of shot noisejHj the 
equivalence of different approaches for calculating noise 
in mesoscopic conductors is not evident. In the regiruepf 
hot electrons the noise suppression factor was founcOS 
to be \/3/4. A gain, this result, which is based on a 
Boltzmann-Langevin equation approach, is restricted to 
quasi-onedimensional conductors. The generalization of 
these results to the case of arbitrary multiterminal con- 
ductors is not obvious. 

We present here the systematic study of transport and 
noise in multiterminal diffusive conductors. This prob- 
lem has been recently addressed by Blanter and Biittiker 
in Ref. |2^, where they use the scattering matrix for- 
mulation followed by an impurity averaging procedure. 
Having the advantage of including quantum phase coher- 
ence, this approach is somewhat cumbersome to general- 
ize to an arbitrary geometry and arbitrary disorder. In 
contrast to this, our approach is based on semiclassical 
Boltzmann-Langevin equation, which greatly simplifies 
the calculations. 

We consider a multiterminal mesoscopic diffusive con- 
ductor (see Fig. |^a) connected to an arbitrary number N 
of perfect metallic reservoirs at the contact surfaces L„, 
n = 1, . . . ,N, where the voltages Vn or outgoing currents 
/„ are measured. The reservoirs are maintained at equi- 
librium and have in general different lattice temperatures 
T„. Unless specified otherwise the conductor has an ar- 
bitrary 3D or 2D geometry with an arbitrary disorder 
distribution. Our goal is to calculate the multiterminal 
spectral densities of current fluctuations at zero 

frequency, u; = 0, 



S'n,n= / dt{SUt)SI„M), 



(1.1) 



where the brackets (. . .) indicate an ensemble average. 
We consider the effects of purely elastic scattering and 
those of energy relaxation due to electron-electron and 
electron-phonon scattering on the same basis. 

Starting our analysis with a brief summary of. the 
Boltzmann-Langevin kinetic equation approachJHjo we 
then apply the standard diffusion app roximation and re- 



1.1 to the solution of a 



duce the problem of evaluating Eq. 
diffusion equation. First, we solve the diffusion equation 
for the distribution function to obtain the multiterminal 
conductance matrix and energy transport coefficicBjts in 
terms of well defined "characteristic potentials" .l3 We 
formulate the Wiedemann-Franz law for the case of a 
multiterminal conductor. Then we turn to the calcula- 
tion of the noise spectrum. We derive the exact general 
formula (3.16) for the multiterminal s pectral density of 
the noise, which together with Eqs. (4^,5.2,3.2 6.5) is 
the central result of our paper. Using this formula we 
demonstrate that the shot- noise suppression factor of 1/3 
is universal also in the semiclassical Boltzmann-Langevin 



approach, in the sense that it holds for a multiterminal 
diffusive conductor of arbitrary shape, electron spectrum 
and disorder distribution. We first prove this for cold 
electrons and then for the case of hot electrons where the 
suppression factor is \/3/4. Thprpby we extend previ- 
ous semiclassical investigationscJo for two-terminal con- 
ductors to an arbitrary multiterminal geometry. This 
allows us then to compare our semiclassical approach 
with the scalieiing matrix approach for multiterminal 
conductors ,lIOE3 in particular with some explicit results 
recently obtained for diffusive conductors.^ The univer- 
sality of shot uoisc proven here gives further support to 
the suggestionnJ that phase coherence is not essentiaLipr 
the suppression of shot noise in diffusive conductors .123 

Another remarkable property of shot noise in meso- 
scopic conductors is the exchange effect introduced by 
BiittikerEZl Although this effect is generally believed to 
be phase-sensitive, we will show that this need not be 
so. Indeed, for the particular case of an H-shaped con- 
ductor (see Fig. |^) we show that exchange effects can 
be of the same order as the shot noise itself even in the 
framework of the semiclassical Boltzmann approach. We 
prove that while the exchange effect measured in differ- 
ent contacts (cross-correlations) can change the sign, it 
is always negative when measured in the same contact 
(auto-correlations). Thus, the auto-correlations are al- 
ways suppressed, in agreement with the Pauli principle. 
Formally, these exchange effects are shown to come from 
a non-linear dependence on the local distribution func- 
tion. Similarly we show that the same non-linearities are 
responsible for non-local effects such as the suppression 
of shot noise by open leads even at zero electron temper- 
ature. 

Finally, we discuss a new phenomenon, namely the cur- 
rent noise in multiterminal diffusive conductors induced 
by thermal transport. We consider the cases of hot and 
cold electrons and prove the universality of noise in the 
presence of thermal transport. We also propose a possi- 
ble experiment which would allow one to measure locally 
the effective noise temperature. Throughout the paper 
we illustrate the general formalism introduced here by 
concrete numbers for various conductor shapes that are 
of direct experimental interest. We note that some of 
results of present paper has been published in Ref. ^ 
in less general form. Here we present the details of the 
derivation of these results and generalize them to a finite 
temperature and an arbitrary electron spectrum (band 
structure). 



II. BOLTZMANN-LANGEVIN EQUATION: 
DIFFUSIVE REGIME 

To calculate the spectral density of current fluctuations, 
we use the Boltzmann-Langevin kinetic equatiorcila 
for the fluctuating distribution function F{p,r,t) = 
/(p,r) + (5/(p,r, i), which depends on the momentum 
p, position r, and time t, 
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eE-dp)F- 



I[F] - I,„^[F] = 6F% (2.1) 



where E(r, t) = E(r) + (5E(r, t) is the fluctuating electric 
field, V = VpE is the velocity of the electron and e is 
its kinetic energy. I[F] = Iee[F] + Ie-ph[F] contains the 
electron-electron and electron-phonon collision integrals, 
respectively (we do not need to specify them here), and 
lim [F] is the impurity collision integral, 

Iim[F] = y ^ {Jp'p ^ Jpp') I 
P' 

Jpp'{v,t) = W^pp.(r)F(p,r,t)[l - i^(p',r,t)], (2.2) 

where the elastic scattering rate from p into p', Wpp'(r), 
depends on the position r in the case of disorder consid- 
ered here. 

The Langevin source of fluctuations SF^{p,r,t) is in- 
duced by the random (stochastic) process of the elec- 
tron scattering which is also responsible for the momen- 
tum relaxation of the electron gas. On the other hand, 
electron-electron scattering conserves total momentum 
of the electron gas and therefore does not contribute to 
6F'^. Furthermore, we neglect the momentum relaxation 
due to electron-phonon scattering and electron-electron 
Umklapp process, assuming that they are weak compared 
to the scattering by impurities in diffusive conductors 
(phonon induced shot noise in ballistic wires has been 
studied in Ref. |3^ ). In other words, we assume that the 
collision integrals /eei^"] and Ie~ph[F] describe only en- 
ergy relaxation process in the electron gas, but it is only 
impurity scattering which gives rise to momentum relax- 
ation and to the shot noise in diffusive conductors. 

To describe the fluctuations SF'^ we make use of 
the Langevin formulation introduced by Kogan and 
Shul'man (Ref. p5| ). In this approach there are two con- 
tributions to the fluctuations of the impurity collision 
integral. First, there is the contribution Iim[Sf] due to 
the fluctuations of the distribution function, which has 
already been included in Eq. (2.2). The second contri- 
bution. Slim [f] , stems from the random character of the 
electron scattering, which is the extra source of fluctua- 
tion 6F'^ occurring on the rhs of Eq. (pj|), i.e.. 



6F' 



p' 



SJpp') , 



(2.3) 



where the random variables SJpp> are intrinsic fluctua- 
tions of the incoming and outgoing fluxes Jpp' . 

Assuming now that the flow of electrons, say, from 
state p to-state p' is described by a Poisson process we 
can writeES 

{SJpp'{r, i)(5Jp^p'^(ri,ii)) 

= S{t - t,)S{r - ri)5pp, Vpi {Jpp'ir, t)) , (2.4) 

where 



Using the preceding two equations together with Eq. 
(^.3|) , we obtain the correlator of the Langevin sources, 

(<5F^(p, r, t)6F%p',r', t')) = <5(i - t')6ir - r') 

X E('^PP' - V'pOW^pp" [/(I - /") + /"(I - ./)] ■ (2.6) 
p" 



with /" = /(p",r), and Wpp-- = I^p, 

Next, we consider the Ihs of Eq. (2T). Since we are 
only interested in the lu — limit of the spectral density 
(the effect of screening on frequency dependent shot noise 
in quasi-one dimensional diffusive conductors has been 
studied recently in Refs. ^ and we may drop the 
flrst term dF/dt in Eq. (2.1). The term eE-dpF can be 
rewritten as follows: eFj-OppF + ev-EdeF, where p_f is the 
momentum at the Fermi surface. From this we see that 
the electric field E induced by an applied voltage plays 
a twofold role: it effects the trajectories and changes the 
energy of electrons. The first effect, eE-^p^F ~ eE/pp, 
is weak compared to v- drF ~ vp/L [L is the size of 
the conductor) and gives contribution of order eV/epj 
which can be neglected.E3 The second effect can be taken 
into account by the replacement e — *■ e — eV{v^t) in the 
argument of the distribution function _F, so that e now 
is the total (kinetic -I- potential) energy of t he e lectron. 
Then, the two terms \-drF + eE,-dpF in Eq. (2T) can be 
replaced by the total derivative v-VF. 

In a next step we apply the standard diffusion approx- 
imation to the kinetic equationEj where the distribution 
function is split into two parts, 

F(p,r,i) =Fo(e,r,i)+l(PF,r)-Fi(e,r,t), (2.7) 

where the vector 1 obeys the equation, 

E M/pp'(r)[l(pf , r) - l(p'^, r)] = v . (2.8) 
p' 

The choice of the distribution function F in the form 
( ^.7|) is dictated by the fact that the impurity collision 
integral Iim[F] does not affect the energy dependence 
of t he d istribution function. Inserting this ansatz into 
Eq. (2.1 ) and averaging subsequently over the momentum 
first weighted with one and then with 1, we arrive at 



V-DFi - I[F] = 0, 



D(VFo + Fi) = \5F^ 



(2.9) 



(2.10) 



Here the overbar means averaging over pi? at the 
Fermi surface inside the Brillouin zone, (...) — 
J dppVp^ (••■)/ / dppVp^, and we introduced the diffu- 
sion tensor, 



D(r) = Dapir) = vJp{pF,r). 



(2.11) 



(Jpp, (r, t)) = I^pp, (r)/(p, r)[l - /(p', r)] . 



(2.5) 



We also used 6F'' — 0, which follows from Eq. (2.6) and 
which reflects the conservation of the number of electrons 
in the scattering process. 
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Using the distribution function ( |2.7| ) we can calculate 
the current density j + Sj = evpD J deFi and due to 
charge neutrality (neglecting accumulation of charge) we 
get the potential, eV + e6V = deFo, where Sc is a 
constant energy near the Fermi level and chosen so that 
— 1, and = J dppVp^ is the density of states 



at th e Fermi level. Upon integration of Eqs. (2^) and 
(2.10) over the energy e the collision integrals vanish and 
we arrive at the diffusion equations for the potential and 
density of current, resp.. 



v-avv = 0, 



-aw. 



(2.13) 



where the conductivity tensor o'(r) = e^h'FD{r) depends 
in general on the position r, and = ei'F J delSF'^ is 
the Langevin source of fluctuati ons of the current density. 
After integrating over e in Eq. (2.6) and averaging over p 



(at the Fermi surface) we use then Eqs. (|2j) and ( |2.1l| ) 



to obtain the correlation function of the Langevin sources 
{Sf^{r,t)Sfp{r'X)) = d{t mr - r')a„^(r)n(r), 



n(r) 



de/o(e,r)[l-/o(e,r)], (2.14) 



where /o is symmetric part of the average distribution 
function / = /o + 1-fi. 



The physical interpretation of Eq. (2.14) is now trans- 
parent: the function 11 describes the local broadening 
of the distribution function and can be thought of as 
an effective (noise) temperature. Then we see that the 



correlator ( 2.14 ) takes an equilibrium- like form of the 
fluctuation-dissipation theorem. This is a direct conse- 
quence of our diffusion approximation. In the diffusive 
regime all microscopic details of the transport and fluc- 
tuation mechanisms are hidden in the same conductivity 
matrix, which appear in the correlator of the fluctuation 
sources (2.14) as well as in the diffusion equations ( ^.12 , 
2.13). It is this fact which leads to the universality of shot 



noise that is independent of microscopic mechanisms of 

the noise. 

Ne xt, su btracting the fluctuating part from Eqs. (2.£) 



and ( 2.10 ) we get the equations for the average distribu- 
tion function /, 

V-aV/o + e'i^F/m =0, /-/o-l-V/o, (2.15) 

which complete the set of coupled equations to be solved. 
Now we speci f y the boun dary conditions to be imposed 
on Eqs. ( [IT^ , 



2.13, 2.15| ). First, we assume that for a 



given energy there is no current through the surface S 
(see Fig. ^). Second, since the contacts with 
perfect conductors the average potential V and its fluc- 
tuations SV are independent of position r on L„ . Third, 
the contacts are assumed to be in thermal equilibrium 
with outside r eservo irs HJ Then we write the boundary 
conditions for ( ^.12 ) and ( 2.13 ), respectively, as 



ds-j(r)|5 = 0, V(r)|i„=K, 



(2.16) 



ds-Sj{r,t)\s^O, SV{r,t)\L„^SVn{t), (2.17) 

and for ( |2T5| ), 

/o(e,r)U„=/T„(£-eK), ds.a(r)V/o(e, r)|s -0, 

(2.18) 

where /T„(e) = [1 + exp(e/r„)] ^ is the equilibrium dis- 
(2.12) tribution function at temperature r„, and ds is a vector 





2.121, 


2.13 


2.17 


, 2.18 


) 



In principle, these equations can be solved exactly which 
would allow us to evaluate for an arbitrary multi- 
terminal geometry of the conductor and for an arbitrary 
disorder distribution. 



III. SOLUTION OF THE DIFFUSION 
EQUATIONS 

A. Multiterminal conductance matrix 

The multiterminal conductance matrix is deflned as 
follows: /„ — "l^m ^nmym (throughout the paper the 
sum over the contacts m runs from m = 1 to m = A^, and 
we omit the limits for convenience). To calculate Gnm 
we ne ed to solve Eqs. ( |2.|ljj ) with boundary conditions 
( ^.16|) . Following Biittikeitf3 we introduce characteristic 
potentials <j>n{^), n = 1, . . . , N, associated with the corre- 
sponding contacts. These functions satisfy the diffusion 
equation and the boundary conditions: 



w-aw<t>n = 0, 



(is-(7V0„ I c = 0, (/)„|l„ 



(3.1) 



(3.2) 



so that they are always positive (/)„(r) > 0, n = 1. . . . , A 
and obey the sum rule (see Appendix ^ , 



(3.3) 



The potential V can be expressed in terms of charac- 
teristic potentials 



(3.4) 



to satisfy the diffusion equation ( 2.12| ) and boundary con- 
ditions (2.16). Then the outgoing current through the 



TO-th contact is Im = ds-j = —J^uIl ds-aS/(j)nVn, 
and using the definition of the conductance matrix we 
get 



4 



Grn.-ri 



(3.5) 



We note here that the multiphcation of the integrand by 
(j}m does not change the integral in the rhs of this equa- 
tion. Moreover, the boundary conditions (B.S) for the 



characteristic potentials allows us to extend the integral 
to the entire surface. Doing so and taking into account 



Eq. (3.1), we then replace the surface integral by an inte- 
gral over the volume of the conductor. We are then left 
with another useful formula for Gnm , 



Gri 



I 



(3.6) 



From this expression and from the sum rule for 0„ it 
immediately follows that Gnm = Gmn, J2rL^nm = 0, 
and Gnn < 0, as it should be. In Appendix ^ we use a 
similar procedure to prove another quite natural property 
of the conductance matrix: Gnm > for n ^ m. 



B. Energy transport coefficients 

We have already seen that the local source of noise 
i s defi ned by the effective noise temperature 11 (see Eq. 
(2.14)), which describes the broadening of the distribu- 
tion function. Another important quantity is given by 
the energy density T(r) acquired by the electron gas due 
to the broadening of the distribution function (effective 
heat density). It is given explicitly by the integral 



UF I dee [fo - 9 {e ~ eV)] = A - ^lypieV)^, (3.7) 



where 9 {e — eV) is the local equilibrium distribu- 
tion function at zero temperature, and A(r) = 
Up deefo{e,r) — h'Fel/2 is the total energy density 

(up to irrelevant constant). 

To calculate T we integrate the first of Eqs. ( 2.15[) over 
e with the weight of e and use the expression ( |3.7] ) for 
A. Then the electron-electron collision integral vanishes, 
and we arrive at the following equation,L3 



V- DVA = V- D\IT + j-E = g. 



(3.8) 



where we introduced the rate of energy relaxation (or 
absorption) due to phonons, ^(r) = —vp deele-ph[f]- 
Eq. (^^) expresses energy conservation: the work done 
on the system by the electric field, j-E, is equal to the 
energy flux to the lattice, g, plus the heati ng of the elec- 
tron gas, —V- _DVT. Integration of Eqs. (2.18) gives us 
the boundary conditions for A, 



A|l„ = Kn = l'F 



ds-DVh\s = 0. 



We assume now that electron-phonon interaction is 
weak (the general case is discussed in Sec. [II C ). Then 
the energy exchange between the electron gas and the 
lattice occurs in the metallic reservoirs far away from the 
conductor, and inside the conductor we hav e q = 0. Eq. 
( |3.8|) for A with the boundary conditions (|3.9| ) can be 
solved in terms of (/>„: A(r)== (/)„(r)A„. Substituting 
this expression into Eq. (3/7) and using Eq. ( ^ ) for V, 
we obtain T, 



T^ufY,' 



Vmf 



(3.10) 



On the other hand, in perfect metallic reservoirs 
(where a oo) the term j-E ~ j^/cr can be neglected in 
Eq. (^.8[). Integration of this equation over the volume 
of the nth metallic reservoir gives the total amount of 
energy transferred to (or absorbed from) the lattice in 
this reservoir, Q„ = J drq{v) = — ds-D\/T. In the 
particular case of thermal equilibrium between the reser- 
voirs, i.e., T„ = r, n = 1, . . . , A'', we can use Eq. ( 3.10| ) 
to get the Joule heat in the nth reservoir. 



2 ^ 



(3.11) 



For a two-terminal conductor, (Vi ~ ¥2)'^ = V'^, G12 = 
G21 = G, we have Qi ^ Q2 — GV'^/2, while the total 
Joule heat is Qi + Q2 = IV. We see in this case that 
the heat contributions releasednOn each side of the two- 
terminal conductor are equaLtJ This general conclusion 
holds for an arbitrary shape of the conductor and arbi- 
trary disorder distribution. This fact is a consequence of 
electron-hole symmetry. 

The following simple analysis of the Eq. ( |3.1l|) exhibits 
its physical meaning. On one hand, the total amount 

of Joule heat, |EmnG'„m(K - V,n)^ = -J2n^nVn = 

Jdrj-E, is simply equal to the total work done by the 
electric field on the system. On the other hand, the value 
^e'^i^FiVn — Vm)'^ can be thought of as the gauge invariant 
difference of energy densities A (i. e. minus the density 
of the potential energy, e'^vpiVm — Vn)Vn) appfied to the 
contacts of the conductor. Then the energy transport 
coefficients, Gnm/^^i^F, are determined by the conduc- 
tance matrix. The last fact is a manifestation of the 
Wiedemann-Franz law, which holds for diffusive conduc- 



tors (together with Eqs. (3.1C) and (3.11)) in the cases 
of cold and hot electrons, as soon as the electron-phonon 
interaction is weak enough. To show the Wiedemann- 
Franz law in its usual form, we consider the thermal 
transport in multiterminal conductors in the absence of 
charge transport, Vn — , n ~ 1,...,A^. In this case 
we can use again the Eq. (3.1C) to calculate the thermal 
current Qn, 



(3.9) 



gg2 



^G 



(3.12) 
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In particular, close to thermal equilibrium, Tm = T 
ATm , we have 



3e2 



^G„™AT„, AT„«T, (3.13) 



where ^^^^Gnm is the thermal conductance matrix. This 
is now the Wiedemann-Franz law in its usual form. 



C. Multiterminal spectral density of noise 

In this section we derive the general formula for the 
multiterminal spectral density of shot noise in the case 
of arbitrary elec tron- phonon interaction. We multiply 
the first of Eqs. (2.13) by V0„ and integrate it over the 
volume of the conductor. Then we evaluate the first term 
in Ihs of the eq uatio n integrating by parts and using the 
second of Eqs. (|2.13| ), / drV(/)„<5j — <f dsd j(l)n- Taking into 



account the boundary conditions (2.17) for Sj and (3.2) 
for we get J drV0„-^j = Sin. Integration by parts in 
the second term of the Ihs of this equation gives / (irV(/)„- 
a'^SV = £ds ■ a\7 (j)nSV = —J2k^rikSVk{t), whe re we 
used Eqs. (|3. 1|, pT ^) for 0„, the boundary condition ( |2.17 ) 



for 6V, 



for the conductance matrix G„ 



This 



leads us to the solution of the Langevin equation ( 2.15 ) 
in terms of characteristic potentials: 

5i„ = (5/„ - Gnr„<5K„ = y drVK -^f- (3.14) 



Now, using the correlator (2.14) for the Langevin sources 
(5j*, we express the generalized multiterminal spectral 
density Snm defined as 



(3.15) 



Sn-m = J dt{6ln{t)6lmm 
— oo 

in terms of characteristic potentials, 



(3.16) 



with the properties: Snra = Smn, S„m ^ 0, and 
Snn > 0. In equilibrium n(r) = 2T, and (3.16) together 
with Eq. lead to the result for the thermal noise. 



S„ 



(3.17) 



which is again a manifestation of the fluctuation- 
dissipation theor em. 

The formula ( 3.16| ) is one of the central results of 
the paper. It is valid for elastic and inelastic scatter- 
ings and for an arbitrary multiterminal diffusive con- 
ductor. The relation of Snm to the measured noise is 
now as follows. If, say, the voltages are fixed, then 
SIn{t) = dln{t), and the matrix Snin — 5*^™ is directly 
measured. On the other hand, when currents are fixed. 



Snm can be obtained from the measured voltage corre- 
lator S^^ = J^^dt{SVn{t)SVmm by tracing it with 
conductance matrices: Snm = Yln 'm' Gnn'Gmm' S^,^,. 
The physical interpretation of ( |3.16| ) becomes now trans- 
parent: n describes the broadening of the distribution 
function (effective temperature) that is induced by the 
voltage applied to the conductor and ctII can thus be 
thought of as a local noise source (see the discussion fol- 
lowing the Eq. (2.14)), while 0„ can be thought of as the 
probe of this local noise. In particular, this means that 
only Snm is of physical relevance but not the current or 
voltage correlators themselves. 

L et us consider now one important application of 
Eq.( |3.16 ). In an experiment one can measure the local 
broadening 11 of the nonequilibrium distribution function 
/o (effective noise temperature 11) at some point r = rp 
on the surface of the conductor by measuring the volt- 
age fluctuations in a noninvasive voltage probe. This is 
an open contact with a small area on the surface of the 
conductor around the point r = Tq. The contact is not at- 
tached to the reservoir so that it does not cause the equi- 
libration of the electron gas , and as a result 11 = const 
around the point Fq. Then, ( 3.16|) can be rewritten as fol- 
lows: 5 = / cirV(/)-(7V(/)n = n(ro) / drV(/)-(TV(/), where (j) 
is the characteristic potential corresponding to the non- 
invasive probe. Using Eqs. ( ^ ) we get S = i?~^n(ro), 
where R is the resistance of the contact which comes 
from the volume around Tq. Finally, taking into account 
( |3.14| , p~l] ) and the fact, that there is no current through 
the voltage probe, SI = 0, we obtain 

S" = i?n(ro). (3.18) 

This means that 11 can be directly measured which 
gives an important information about nonequilibrium 
processes in the conductor. Eq. (3.18) resembles the 
fluctuation-dissipation theorem. This is so because there 
is no transport through the noninvasive probe, and there- 
fore one can think of the probe as being in local equilib- 
rium with the effective temperature 11. For this reason 
our consideration restricted to the diffusion regime can 
in principle be applied to the case of the tunnel coupling 
between the probe and conductor. A possible experiment 
that could measure shot noise at local tunneling contacts 
is discussed in detail in Ref. The above result can 
be easily generalized to take into account the equilibra- 
tion by the contact (see Sec. Vl). There w ill be then an 



additional noise supp ression factor in Eq. (3TS 



2.1C 



2.1J 



We note that (|3.16| ) together with Eqs. (|2^ 
for the average distribution function / and Eqs. (3.1 
3.2 ) for the characteristic potentials can serve as a start 
ing point for numerical evaluations of Snm- For purely 
elastic scattering as well as for hot electrons it is even 
possible to get closed analytical expressions for Snm as 
we will show next. The physical conditions for different 
transport regimes are discussed in Ref. |l|. In Sec. IV 
and [V| we will consider the charge transport (T„ = T, 
n = 1 . . . N), and in the Sec. VI we will discuss the ther- 
mal transport {Vn =0, n = 1 . . .N). 
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IV. ELASTIC SCATTERING 

In the case of purely elastic scattering, /[/] = 0, the 
average distribution function satisfies the diffusion equa- 
tion 



V-<7V/o = 0, 



(4.1) 



and the boundary conditions ( 2.18| ) with T„ = T (i. e. 
in the charge transport regime) . Using this equation one 
can prove (see Appendix ^ that for elastic scattering 
cross correlations (n ^ m) are always negative, in agree- 
ment with the general conclusion of Ref. 

Eq. (4.1) can be solved in terms of (/)„: /o = 
(jj ufris — eVn). Substitutin g th is solution into Eq. 
(2.14) and using the sum rulp_.(^3|) for we can ex- 
press n in the following form, til 

n = 2 f deJ^^k^iMe - eVk)[l - Me - eVi)]. (4.2) 

•' k,l 



Performing the integration over e we obtain, 

'e{Vk-Vi) 



n = e ^ (pkipi (Vfc - Vi) coth 



k.l 



2T 



(4.3) 



which in combination with Eq. ( ^.16| ) gives the final ex- 
pression for Snm which is valid for purely elastic scatter- 
ing. Eq. (4.3) describes the crossover from the shot noise 
in multiterminal diffusive conductors (T — > 0), 



Il = eJ2^kMVk-Vi\, 



(4.4) 



k.l 



to the equilibrium Johnson-Nyquist noise given by ( 3.17 ) 



A. Universality of noise 

Now we are in the position to generalize the ffifj^ of 
universality of the 1/3-suppression of shot noiseLTllB to 
the case of an arbitrary multiterminal diffusive conduc- 
tor. To be specific we choose Vn = 0, for n 1, i.e. only 
contact n — 1 has a non- vanishing voltage. Then, using 
the sum rule (3.3) for (/)„ , we get 



n = 2e0i(l - 0i)FiCoth(eVi/2T) 
-|-2T(1 - (Pif + 2T(f>l . 



To get Sin we substitute this equation into ( 3.16 ) and 
evaluate the first term as follows: / drV(/)„-tTV(/)i(/)i(l — 
(pi) — f ds-a\/(j)n{(j)l/2 — (j)l/3) = — Gi„/6, where we used 
( |3.1| , 3.2). Similarly, for the integrals in the second and 
third term we get: J (irV0„-(TV(/>i(l — = J dx'^<l>n- 
a"^ 4'i4>\ — — Gi„/3. Combining these results we arrive 
at 



^i„ = -^Gi„[4r 



sViCoth {eVi/2T)] 



(4.5) 



When Vi = we get S'i„ = — 2Gi„T, and the formula for 
the Johnson-Nyquist noise is recovered. When T = 0, we 
express 5'i„ in terms of outgoing currents, /„ = GinVi: 



Sin 

Sii 



1 



-elln 



n^l. 



(4.6) 



We note that above derivation is valid for arbitrary im- 
purity distribution and shape of the conductor, and for 
an arbitrary electron spectrum (band structure). In 
this sense the suppression factor i is indeed universal. 
This generalizes the known universality of a two-terminal 
conductoill3 to a multiterminal geometry. 

Finally, we mention here some inequalities (derived in 
Appendix ^ , which can be used to estimate the spectral 
density Snm in the T = limit. First, the correlations 
are bounded from below, 



Snn — g^l^nli 



(4.7) 



but due to th e nonlocality of the noise (see the discussion 
in Sec. iVC) there can be no upper bound in terms of 
the current /„ through the same contact. In other words, 
the current /„ flowing through the n-th contact creates 
the noise ^e/„ in this contact. However, other contacts 
also contribute to the noise in the n-th contact, and this 
contribution is not universal and makes the noise arbi- 
trarily larger compared to the value |e/„. Nevertheless, 

kj: 



we can write: 11 < max{|V/j - 
gives the rough estimate 



Vi\}, 



1, 



Snn < e|G„„| max{|Vfc - Vi\}. 



, iV, which 



(4.8) 



In contrast, the cross correlations possess an upper 
bound, 



|5'nm| — ^{Snn ^ Snini)- 



(4.9) 



Snm vanishes when the n-th and m-th contacts are com- 
pletely disconnected. 



B. Wide and star-shaped conductors 

Next we specialize to two experimentally important 
cases. First we consider a multiterminal conductor of 
a star geometry with N long leads (but with otherwise 
arbitrary shape) which join each other at a small cross- 
ing region (see Fig. |l]b). The resistance of this region 
is assumed to be much smaller than the resistance of 
the leads. In the second case the contacts are connected 
through a wide region (see Fig. |l]c), where again the re- 
sistance of the conductor comes mainly from the regions 
near the contacts, while the resistance of the wide region 
is negligible. 
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Both shapes are characterized by the requirement that 
w/L <^ 1, where w and L are the characteristic sizes of 
the contact and of the entire conductor, resp. In both 
cases the conductor can be divided (more or less arbi- 
trary) into N subsections Tk, k = 1,...,N, associated 
with a particular contact so that the potential V is ap- 
proximately constant (for w/ L <^ 1) on the dividing sur- 
faces Ck ■ Each subsection then can be thought as a two- 
terminal conductor with the corresponding characteristic 
potential Ok{r), 

y-aVOk^O, ds-aV9k\s^ek\L,^0, ^fclc^l- (4.10) 

We will show now that both, the multiterminal conduc- 
tance matrices Gnm and the spectral densities Snm, can 
be expressed in terms of the conductances Gk of these 
subsections, 

Gk = - j ds-aVOk = j ds-aVOk- (4.11) 

Since each potential (/>„ is approximately constant in 
the central region of the multiterminal conductor, we can 
write, 



'^n(r)|cfc = a„ = consi., ^a„ = l. 



(4.12) 



for an arbitrary fc = 1, . . . , A^, where the second equation 
follows from the sum rule for (^„. Com paring Eqs. (3.1, 



3.2, 4.12 ) with the definition of 9k ( 4.10 ), we immediately 
obtain 



a„6lfe(r) + [l-e'fc(r)]5„fc. 



(4.13) 



The calculation of Gnm and Snm is now straightfor- 
ward. We substitute ( f4.13 ) into (^) and use Eq. ( |4.1l| ) 
to get 



Gn 



)G 

n 7 



Olfn — Gin/ 



(4.14) 



where G = G„, and the equation for a„i follows from 
X^ri ^rmi = 0. Substituting ( 4.13| ) into ( 3.16| ) and apply- 
ing similar arguments as above in the proof of the 1/3- 
suppression we find the explicit expressions (for details 
of the derivation see Appendix ^) 

Snm — ^ ^ ^n^k{,Jk ~^ Jn){,^nm ^k7ti) '^Gjijn'-^j 



'^» = IlGz(K-V,)coth 



eiVn-Vi) 



2T 



(4.15) 



We note t hat t his result is a consequence of above approx- 
imation ( [4.12| ). Comparing the resistance of the subsec- 
tions to the resistance of the central region of the con- 
ductor (which is neglected) we find that the corrections 
to Eq. (4.12) and consequently to Eq. ( 4.15| ) are of or- 
der w/L in 3D and for a star geometry in 2D, and up to 



corrections of order [ln(L/w)] for wide conductors in 
2D. 

In principle, ( 4.15| ) and (4.14) allow us to calculate 
the noise for arbitrary voltages and temperature, but 
for illustrative purposes we consider the simple case of 
a cross-shaped conductor with four equivalent leads (i.e. 
a„ — 1/4) and T = 0. Suppose the voltage is ap- 
plied to only one contact, say Vi > 0, V n^i = 0, and 
I ^ -h ^ 3/„^i > 0. Then, from ( |4T5| ) we ob- 
tain: S'li = |e/, = S'ls = 5*14 = 
being in agreement with the universal 1/3-suppression 
proven above. Then, S22 — = S44 = |e/, and 
S23 ~ M24 — S34 = —jgcl. These numbers seem to 
be ncwcS and it would be interesting to test them exper- 
imentally. 



C. Nonlocality and exchange effect 

We are now in the position to address the issue of non- 
locality and exchange effect in shot noise {T — 0) in mul- 
titerminal conductors. For this we consider for instance 
a star geometry and assume that the current enters the 
conductor through the n-th contact, i.e. /„ = — J, and 
leaves it through the m-th contact, i.e. Im = I, while the 
other contacts are open, i.e. Ik = for k ^ n,m. From 
( [4.14 ) we obtain for the conductance GnGm/ {Gn + Gm) 
(two contacts are in series), and we see that it does not 
depend on the other leads, which simply reflects the lo- 
cal nature of diffusive transport. However, contrary to 
one's first expectation, this locality does not carry over 
to the noise in general. Indeed, from ( 4.15| ) it follows 
that Snm = —\{oi-n + ctm)^/. The additional suppres- 
sion factor < q;„ -|- am < 1 for iV > 2 reflects the 
non-locality of the current noise. For instance, for a cross 
with = 4 equivalent leads we have a„i = a„ = 1/4, 
and thus Snm. = ^^^J. An analogous reduction factor 
was obtained in Ref. ^ under a different point of view. 
Hence, one cannot disregard open contacts simply be- 
cause no current is flowing through them; on the con- 
trary, these open contacts which are still connected to 
the reservoir induce equilibration of the electron gas and 
thereby reduce its current noise. We emphasize that this 
non-locality is a classical effect in the sense that no quan- 
tum phase interference is involved (phase coherent effects 
are not contained in our Boltzmann approach). On the 
other hand, the origin of this non-locality can be traced 
back to the non- linear dependence of H on the distri- 
bution / in ( 2.14 ), which is a consequence of the Pauli 
exclusion principle. _ 

Next we discuss exchange efFect£0 in a fouc, terminal 
conductor. According to Blanter and Biittikerc^l they can 
be probed by measuring S13 in three ways: Vn — Vo(5n2 
(A), Vn = VoSn4 (B), and K = Vo5n2 + VoSn4 (C). Then 
we take AS'13 = <5'f^ — S'j^ — 5^ as a measure of the 
exchange effect. This experiment is analogous to the^x- 
periment of Hanbury Brown and Twiss in opticso It 
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measures the interference of electrons coming from mu- 
tually incoherent sources, which is caused by the indistin- 
guishability of the electrons. Naively, one might expect 
that this interference effect averages to zero in diffusive 
conductors. However, it comes now as some surprise that 
in our scmiclassical Boltzmann approach AS'13 turns out 
to be non-zero in general and can even be of the order 
of the shot noise itself. Again, the reason for that is 



that n is non-linear in /o (see (2.14)), which is the con- 
sequence of the Pauli exclusion principle. So, the value 
n'-^ — 11'^ — which enters AS'13 is not necessarily zero. 
Indeed, while exchange effects vanish for cross shaped 
conductors (in agreement with Ref. ^ up to corrections 
of order w /L which are neglected in our approximation) , 
it is not so for an H-shaped conductor (see Fig.0d). Cal- 
culations similar to those leading to (4.15) give for this 



AS13 = 



1 eVoG^Go 
24 (G + 4Go) 



(4.16) 



where Gn — G/A are the conductances (all being equal) 
of the outer four leads, while the conductance of the 
connecting wire in the middle is denoted by Gq. This 
exchange term A5'i3 vanishes for Go — > 00, because 
then the case of a simple cross is recovered, and also 
for Go 0, because then the 1-st and 3-rd contacts 
are disconnected. AS'13 takes on its maximum value for 
Go = G„ and becomes equal to ^e/^, where is the 
current through the 2-nd contact for case (A). 

Although AS13 is positive in the example considered 
above this is not the case in general. For an arbitrary 
four-terminal geometry of the conductor the exchange 
effect can be expressed in terms of characteristic poten- 
tials: 



AS„ 



-46^0 / drS/(j)n-a\7(j). 



(4.17) 



where all indices are different. From this general for- 
mula it follows that AS„m = AS^;, and AS„m -I- AS„; + 
AS„fc = 0. The last equation means that the exchange 
effect can change sign, i.e. cross correlations can be either 
suppressed or enhanced. 

On the other hand, the set-up can be slightly modi- 
fied: instead of cross correlations, the noise density in 
one of the contacts of a multiterminal (A'' > 2) diffusive 
conductor is measured, say Sn, while the electrons are 
injected through the contacts 2 (A), 3 (B), and 2 and 3 
(C). Again, ASn = Sfi - Sf\ - Sfi is a me asure of the 
exchange effect. Then, it follows from ( 1.17| ) that 



46^0 / drV0i-(TV0i02(/>3 < 0, (4.18) 



ASn 



i.e. the correlations are always suppressed due to the in- 
terference effect, which is a direct manifestation of the 
Pauli principle. In the particular case of star-shaped con- 
ductors we have 



ASn 



4 . G1G2G3 



G 



N 

E 

n=l 



Gn 



(4.19) 



The suppression of noise due to the interference of mu- 
tually incoherent electrons was recently observedpin an 
experiment with a ballistic electron beam splitter JlJ We 
have shown here that this effect is also observable in 
mesoscopic diffusive conductors. 



V. HOT ELECTRONS 

We consider now the case of "hot" electrons where 
lee 7^ 0, but still le-ph = 0, and we assume that electron- 
electron scattering is sufficiently strong to cause thermal 
equilibration of the electron gas (i.e. lee = VDrle ^ 
where D is the diffusion coefficient and Tee the electron- 
electron relaxation time). The average distribution then 
assumes the Fermi-Dirac form: 



/o(e,r) 



1 + exp 



eV{v) 



Teir) 



(5.1) 



with the loc al ele ctron temperature Te(r). Substituting 
this fo into ( 2.14 ) we imm edia tely get 11 (r) = 2Te(r). On 
the other hand, from Eq. (|^ it follo ws th at (Teir))^ = 
(6/7r2zyp.)T(r), where T(r) is given by ( |3lo| ) with T„ = T 
(i. e. in the charge transport regime). Thus, we finally 
obtain 



U = 2T 



2 0n (I'm {Pn-PmY 



V3eK 



2ttT 



(5.2) 



which in combination with ( |3.16 ) gives the general so- 
lution for the case of hot electrons. We would like to 
note here that the cross correlations are always negative 
also in the case of hot electrons (the proof is given in the 
Appendix 

A. Universality of noise 

Next we show that the shot-noise suppression factor 
\/3/4 (Refs. |l|,||) for hot electrons in a multiterminal 
conductor is also universal. As before we can consider 
the case where the voltage is applied to only one contact: 
14 = ViSni- Then 



n = 2r^l + 4/32 (</,!- <^2)^ 
where (3 = (3i. Using the relation 

2^1 - $2V$ = V |arcsin$ + <i>y/l - ^^j ^ 



(5.3) 



$ = /3 1 + /?^ 



(201 - 1), (5.4) 
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we transform the volume integral in (3.16) into a surface 
integral and obtain the spectral density of noise: 



'5'lii — —GinT 



1 



f3+— ] arctan/? 
P 



/3 = 



27rT 



(5.5) 



This expression describes t he cr ossover from the thermal 
noise (/3 <C 1) given by ( 3.17 ) to the transport noise 
iP » 1) 



Sin — ~~^s\In 



n^l, 



(5.6) 



This general result shows that in the case of hot elec- 
trons the shot-noise suppression factor •\/3/4 is indeed 
universal, i.e. it does not depend on the shape of the 
multiterminaL diffusive conductor nor on its disorder 
distribution.c3 

The origin of this universality becomes clear from the 
following argumentation. We have seen that the distri- 
bution of the effective noise temperature n(r) for the 
case of hot electrons i s co ntrolled by the transport equa- 
tions for the energy, (3^, 3^) through the heat density 
T(r). The spectral density of noise, in turn, is giv en by 
n(r) through the transport equations for charge, (2.13, 
2.17 ). On the other hand, according to the Wiedemann- 
Franz law both, the energy and charge transports, are 
determined by the same kinetic coefficients, namely, by 
<j{y). Thus, the physical origin of the universality of 
the suppression factor -^3/4 can be traced back to the 
Wiedemann-Franz law. Conversely, a violation of the 
Wiedemann-Franz law will cause deviations from univer- 
sality. 

We would like to note here that the universality of the 
noise (for cold and hot electrons) has been proven here 
for the case where the voltage is applied to only one con- 
tact of a multitcrminal diffusive conductor which made it 
possible to express the spectral densities Snm in terms of 
conductances Gnm- This is no longer possible in general. 
Nevertheless, in the case of a 2D geometry and isotropic 
conductivity, Oapir) = a{r)6af3, both G„m and Snm are 
of the same universality class. Indeed, one can easily see 
that they are invariant under conformal transformation 
of coordinates. Gnm and Snm can be expressed in terms 
of characteristic potentials </>„, which satisfy the confor- 
mal inva rian t diffusion equation (3.1) and boundary con- 
ditions (3.2). Moreover, the combination (irV0„-V0„i 
does not change with the conformal transformation of 
coor dinates, which fina lly makes the integrals for Gnm 
(3.6) and for Snm ( 3.16 ) conformal invariant. 

We close this section by another illustrative example. 
Let us consider again a cross-shaped conductor with four 
equivalent leads,L2l Gn = G/4, at T = and where we 
choose Vn = ViSin, I = —Ii = 3/„^i > 0. We then 



find Sii = -^e/, and Sin = 't- 7^ li while 

o _ / 35^/3 _ J2 \ pT and <? ^ - - ( 13\/3 _ J_\ 
Onn — y log 3-K ) '^'"^ 'Jn^m — ^ iqS 'iTT ) ^ 

for n^m ^ 1. These new numbers are consistent with the 
universal factor ^3/4. 



VI. NOISE INDUCED BY THERMAL 
TRANSPORT 

In this section we address a new phenomenon, namely 
the current noise in multitcrminal diffusive conductors in 
the presence of thermal transport. We assume no energy 
relaxation in the conductor due to phonons, le-ph — 0, 
and no voltage is applied to the contacts, Vn = 0, 
n = I, . . . , N, which are kept in equilibrium at different 
temp eratu res r„ . The thermal transport is considered in 
Sec. [II B, where the outgoing thermal currents Q„ are 
calculated (see Eqs. ( p.l2D and (|3.13D). We turn now to 
the calculation of the spectral density of noise. 



A. Elastic scattering 



To calculate 11 we need to know the distribution func- 



tion /o. It o beys the Eq. (4.1) with the boundary con- 
ditions (2.18) in the contacts. The solution then reads 
explicitly, 



/o 



and with the help of 



(6.1) 



we get, 



oc 

Z..=T.T,/rf.[l-tanh(T„s)tanh(T..)], 



(6.2) 



This together with Eq. (|3.16| ) gives the spectral density 
of noise Snm ■ 

In equilibrium r„ = T, Z^i = 2T and (|6.2|) and we find 



for the equilibrium noise, Sn 



-2GnmT. On the other 



hand, if for example Tk'^Ti, then Zu = {2\n2)Tu. We 
consider then two situations, where e.g. either Ti = T 
and T„^i = 0, or Ti = and T„^i = T. In other words, 
only one contact is cither heated up to high enough tem- 
perature T or cooled down to zero temperature and the 
other cantacts are kept at the same temperature. Using 
the sum rule for (/)„ and carrying out the integration in 



(I 



16) we obtain for both cases, 
2 



Sir 



-(l + ln2)Gi„r. 



(6.3) 



Sin for this two situations can be expressed in terms of 
thermal currents Qn 

^i„ = ±4(l + ln2)(e/7r)2T-iQ„, (6.4) 

with the sign depending on the sign of Q„. 
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B. Hot electrons 

We consider now the case of hot electr ons, where 11 — 
2Te, while ^ J2k^kTl (see Eq. (|3lo|)). Substituting 
this into ( 3.16| ) we get, 



Snm = 2 / drS/(j)n-a-\7(f>. 




(6.5) 



In particular, if the electron gas in the conductor is 
pushed out of equilibrium by heating (or cooling) one 
of the contacts (with n ~ 1) while the other contacts are 
kept at the temperature Tn = T2, n ^ 1, the integral can 
be calculated explicitly, and we have 



4 

Sin — —-^Gin — rji 

6 ^1 + ^2 



T1T2 



(6.6) 



In the cases where Ti = T > and T2 = T > Ti , we 
obtain with the help of Eq. ( 3.13| ) that the spectral den- 
sity Sin can be expressed in terms of thermal currents: 



Sm - ±8(e/7r)2r-iQ, 



(6.7) 



Expressions (p^, 6.7) are analogous to Eqs. (4.6, |5.6|) 



and reflect the universality of the noise in the presence 
of thermal transport. 



VII. CONCLUSION 

In conclusion, we have systematically studied the 
transport and noise in multiterminal diffusive con- 
ductors. Applying a diffusion approximation to the 
Boltzmann-Langevin kinetic equation we have derived 
the diffusion equations for the distribution function and 
its fluctuations. We then solved these equations in gen- 
eral terms of well defined "characteristic potentials" and 
we derived exact formulas for the conductance matrix, 
energy transport coefficients and the multiterminal spec- 
tral density of noise. In this way we have obtained the fol- 
lowing results. In both regimes of cold and hot electrons 
the shot noise turns out to be universal in the sense that it 
depends neither on the geometry of a multiterminal con- 
ductor and the spectrum of carriers, nor on the disorder 
distribution. We have studied the noise in the presence of 
thermal transport and find that being expressed in terms 
of thermal currents it is also universal. We believe that 
the origin of this universality lies in the fact that in the 
diffusive regime the correlator of the local current den- 
sities (Langevin sources) takes an equilibrium-like form 
of the fluctuation-dissipation theorem involving an effec- 
tive noise temperature. Thus, the transport and noise 
properties are determined by the same conductivity ten- 
sor. One can surmise then that the proven universality 
holds as long as the energy transport is governed by the 
Wiedemann-Franz law. 



The exchange effect is proven to be non-zero even 
within our semiclassical Boltzmann approach. The ex- 
change effect can change sign when measured in cross- 
correlations, and (in agreement with the Pauli princi- 
ple) it gives always negative contribution to the auto- 
correlations. The exchange effect comes from a non-linear 
dependence on the local distribution function. Similarly 
we show that the same non-linearities are responsible for 
non-local effects such as the suppression of shot noise by 
open leads even at zero electron temperature. 

Finally, we have proposed a possible experiment which 
would allow one to locally measure the effective noise 
temperature, and we have given new suppression factors 
for shot noise in various geometries which can be tested 
experimentally. 
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APPENDIX A: 

In this Appendix we derive some properties of the char- 
acteristic potentials (/)„, multiterminal conductance ma- 
trix Gnm and spectral density of shot noise Snm ■ 

First, we show that </)„ > in the conductor. Accord- 
ing to the boundary condition ( |3.2| ), being negative (/)„ 
would take on its minimum value at some point r = Tq 
in the conductor. If this happened inside the conductor, 
then V(j!)„(ro) = 0, and V • (T(ro)V</)„(ro) > 0, because 
(7 is a positive defini te m atrix. This, however, would 
then contradict Eq. ( ^.l| ). If (/)„ took on its minimum 
value on the open surface of the cond ucto r (ro G S) , then 
V||(/)„(ro) = 0, and according to Eq. ( |3.2| ) V_L(/)„(ro) — 0. 
Aga in, V • (j(ro)V0„(ro) > in contradiction with Eq. 
( |3.lD . Thus, we see that the characteristic potentials can- 
not be negative, 4'n > 0. 

The sum rule (3.3) for 0„ follows from the observation 
that the function 0(r) = 1 is a unique solution of the 



diffusion equation V • aVcj) — with the boundary con- 
ditions ds-(7V(/)|s — 0, and = 1, Vm. It follows 
from Eqs. ( |3.1| , 3.2) that the function obeys 



the same equation and boundary conditions, and there- 
fore (j) = J2n 't'ni'lc) = 1. This sum rule can also be de- 
duced from the fact that physical observables are invari- 
ant under a global shift of the energy scale by a constant 
value .cj 

Now we prove that Gnm > for n ^ m. W e note 
that the multiplication of the integrand in ( |3.5| ) by 0^ 
and extension of the integral to the whole surface does 
not change the integral, so that Gnm = ~ § ds-aW4>n<Pm- 
Then, using the Eq. (^J) we replace the surface integral 
by the integral over the volume of the conductor, Gnm = 
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—2 J dr'V(j>n-d''V(j)m'Pm- Finally, we calculate this integral 
by p arts and take into account the boundary conditions 
( |3.2D for (j)n to get for the conductance matrix, 

Gnm = '2 J drV(j)jn-d-V(t)m4>n, u ^ m, (Al) 

from which it follows that G„m > 0. Interestingly, a 
similar procedure for n ~ m gives. 



Gr. 



2 / (irV(/)„-(TV(/)„ 



(A2) 



and Gnn < 0. 

Now we prove that in the case of elastic scattering 
cross correlations Snm , n ^ m, are alway s negative. We 
calculate Snm in (3.16) with 11 from (2.14) integrating by 
parts, 

Snm = ds-&{\/(l)n<f),n + V (/)„ ) 11 

- i y" (irVn-(T(V(/)„0,„ + V0„(?!)„). 



From ( ^.18 ) it follows that only conta ct s urfaces con- 
tribute to the first integral, and with ( |3.5[ ) we get for 
the first term: —Gnm{Tn + Tm)- In the second term we 



again calculate the integral by parts and use (|2.14| ) to 
get: f d s-(7V n(/)„(/)m-2 // drd£0„0mV/o-CTV/o. Ac- 
cording to (2.18) and ( |3.2| ) the surface integral disappears 
and we arrive at the following result 



Sn 



Gnm{Tn + T„) - 2 / / drd£(/)„0„ V/q • itV/o 



(A3) 



can put T = in the following. Then, one can easily see 
that 

Il = eJ2^kMVk-Vi\>2eq}nJ2(j)i\Vi-Vn\ 



2e^n\V-Vn\. 



We substitute this into ( |3.16[ ) and write another set of 
inequalities, 

Snn >2ej drV?!)„-<7V0„(/)„|y - K| 

>2e J dr\/(j)n-a\7(t)nMV -Vn) 

Integrating by parts, 

2 J drV(/)„-<7V0„0„(y-K) 

= (j) ds-a\/(t)n(j)l{V - Vn) - J dr\/V-a\/(j), 



we obtain expression (4/7). The inequality Eq. ( |4.9| ) 
immediately follows from 

2|V0„-(tV(/),„| < V(/)„-ctV(/)„ + V(/),„-ctV(/),„ , (A5) 

and evidently holds also for inelastic scattering. 

APPENDIX B: 



for n 7^ TO, from which it follows that the cross correla- 
tions are negative. 

We apply similar arguments to prove that the cross 
correlations are always negative also in the case of hot 
electrons. We calculate the integral in Snm — 2 J drW(f>n- 
aVipmTe and use the boundary conditions for 0„ and Tg 
to get 

Snm = -Gnm{Tn + T„i) + y (ir(/)„0,„ V • aVTe . 

Then we use (Te(r))2 = {6/n^iyF)T{r) and Eq. (U) (g = 
0) to write TeV-o-VTe = -WTe-aWTe-3{e/Tr)'^WV-aWV. 
Substituting this into above equation we obtain 

Sn7n — GnmiTn ^" Tjn) 

- J drcbn(l)mT-^[yTe-aVTe + 3{e/iTfVV-aVV], (A4) 

which shows the nega tivity of cross correlations. 

The inequality (4.7) can be derived as follows. First we 
note that according to Eq. (4.3) correlations grow with 
temperature. Therefore, without loss of generality we 



Th e derivation of ( 1.15| ) proce eds a. s follows. We use 
(4.13) to replace the integral in ( 3.16| ) over the volume 
of the conductor by the sum of integrals over subsections 
Ffe, 

Snm ='^{an-Snk){am-Smk) J drV6k-^V6kU. (Bl) 



Using the contraction 

Zig^e{Vi~Vq)coth 
we express 11 in terms of 6k, 



eiVi - V,) 



2T 



(B2) 



n|r, ^J2^aiek + ii-ek)dik] 

Iq 

X [aqOk + (1 - 9k)5qk]Zlq. 



(B3) 



Then we substitute this expression into Eq. (Bl ), ca lcu- 
late the integrals over Ffc with the help of Eqs. (4.11), 
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/ 



drV9k-aV0kOl = 2 J drV0fc-<TV0fe0fe(l - Ok) 



(B4) 



and use the symmetry Ziq = Zqi to get 

Snm = o X/'-'^" ~ ^nfc)(am - (^mfc) 



klq 



X (a/a, + a/(5qfc + 5ik5qk)GkZiq. 



(B5) 



Finally, introducing the notation, Jq = e~^^jGi^g;, 
and using Z^^ = 2T and = Gk/G, we c arry out the 
summation over k in (B5) to arrive at (4.15) for Snm- 
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